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ABSTRACT 

We present the concentration (c)-virial mass (M) relation of 39 galaxy systems ranging in mass 
from individual early-type galaxies up to the most massive galaxy clusters, (0.06-20) x 10 14 M Q . We 
selected for analysis the most relaxed systems possessing the highest quality data currently available 
in the Chandra and XMM public data archives. A power-law model fitted to the X-ray c — M relation 
requires at high significance (6.6a) that c decreases with increasing M , which is a general feature 
of CDM models. The median and scatter of the c — M relation produced by the flat, concordance 
ACDM model (il m = 0.3, eg = 0.9) agrees with the X-ray data provided the sample is comprised 
of the most relaxed, early forming systems, which is consistent with our selection criteria. Holding 
the rest of the cosmological parameters fixed to those in the concordance model the c — M relation 
requires 0.76 < < 1.07 (99% conf.), assuming a 10% upward bias in the concentrations for early 
forming systems. The tilted, low-cg model suggested by a new WMAP analysis is rejected at 99.99% 
confidence, but a model with the same tilt and normalization can be reconciled with the X-ray data by 
increasing the dark energy equation of state parameter toi«~ —0.8. When imposing the additional 
constraint of the tight relation between <jg and Q m from studies of cluster abundances, the X-ray 
c—M relation excludes (> 99% conf.) both open CDM models and flat CDM models with 17 m » 1. 
This result provides novel evidence for a flat, low-f2 m universe with dark energy using observations 
only in the local (z <C 1) universe. Possible systematic errors in the X-ray mass measurements of a 
magnitude » 10% suggested by CDM simulations do not change our conclusions. We discuss other 
sources of systematic error in the measurements and theoretical predictions that need to be addressed 
for future precision cosmological studies using the c — M relation. 

Subject headings: X-rays: galaxies: clusters — X-rays: galaxies — dark matter — cosmological pa- 
rameters — cosmology: observations 



1. INTRODUCTION 

The current cosmological paradigm of a flat uni- 
verse with low matter density (O m ) consisting of Cold 
Dark Matter (CDM) and dark energy has become es- 
tablished recently through a variety of astronomical 
observations, especially the Cosmic Microwave Back- 
ground (CMB), supernova Hubble diagram, and statis- 
tics of large galaxy surveys (e.g.. ISpergel et all 120031 : 
lAstier et all 120061: ISpergel et al.ll2006fh These key tech- 
niques probe the large-scale, high-redshift universe where 
the development of cosmic structure may be studied in 
the linear clustering regime (5 <C 1). Although the stan- 
dard mo del appears to have withstood every test thus 
far (e.g.. iPrimackl l2006h . testing the ACDM model on 
small scales, fully into the non-linear regime (S ^> 1), 
will further our understanding of galaxy formation and 
evolution and may also refine the fundamental parame- 
ters of the world model itself. 

The dark matter halos of galaxies and clusters are the 
largest virialized manifestations of small-scale, non-linear 
clustering. An important probe of the structure of dark 
matter halos is the relationship between concentration 
and virial mass. The concentration parameter is defined 
as, c = r\/r Sl where ta - the "virial radius" - is usually 
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taken to be the radius within which the average density 
equals Ap c , where p c is the critical density of the uni- 
verse, and A is a number typically between 100 — 500. 
The quantity r ? is t he scale radius of the NFW profile 
([Navarro et al.lll997l ). but it is replaced by r_ 2 , the ra- 
dius where the logar ithmic density slope equals -2, for 
more general profiles (|Navarro et al.|[2004l : I Graham et al.l 
2006). The virial mass (M) is the mass enclosed within 

N-body simulations of CDM models find that c declines 
slowly with increasing M bu t with substantial i ntrin- 
sic scatter independent of M (|Bullock et al.ll200ll here- 
after BO l^Dolag et al.ll2004 hereafter D04:lKu"hlen et all 
120051: IShaw et alj 120061: iMaccio et al.ll2006l) . For a given 
M the value of c in CDM models varies significantly 
with changes in cosmological parameters, particularly 
the normalization f2 m , and w, the para meter rep- 
resen t ing the dark ener gy eq uation of state (lEke et al.l 
12001b lAlam et al.l 120021 : D04: iKuhlen et al.ll2005l ). The 
amount of intrinsic scatter in the c — M relation is a 
robust prediction of CDM models, though preferentially 
selecting relaxed, early forming systems yie lds smaller 
scatter than for the entire halo popula tion ( Jing 120001 : 
IWechsler et al.l 120021: iMaccio et al.ll2006l ). 

The expected decrease of c with increasing M in CDM 
models has yet to be confirmed by observations. Op- 
tical studies of groups and clusters (M > 10 13 M Q ) 
using several diffe r ent t e chniques - galaxy dynamics 
(jBiviano fc Saluccil [200 6: Loka s et al l l2006h . redshift- 
space caustics ( Rines fc Dia fcrio 200(f), and weak gravi- 
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tational lensing (|Mandelbaum et"aill2006l ) - all are con- 
sistent with no variation in c. However, the c — M re- 
lations obtained by these studies are consistent with the 
concordance ACDM model. 

X-ray studies of t he c — M relation using data from 
the ASCA satellite flWu fc Xud 1200(1 iSato et alJl200Ch 
found, c ~ -^20o' 5 , m uch steeper than pr oduced in CDM 
models (B01; D04; iMaccio et all 120061) . However, it 
has only become possible with data from the Chandra 
and XMM satellites to obtain spatially resolved temper- 
ature measurements of sufficient quality for reliable X- 
ray constraints on the mass profiles of elliptical galax- 
ies, galaxy groups, and c lusters (e.g., see reviews by 
iBuotd l2004t I Arnaudl 120051 . and references therein). In- 
deed, our recent study of the mass profiles of seven 
early-type galaxies with C handra indicates c — M val- 
ues consistent with ACDM ([Humphrey et al.l l2006) . The 
c — M relatio ns for clusters (M > 10 1 4 M Pl ) obtained 
using XMM ([Pointecouteau et all [2005) and Chandra 
(| Vikhlinin et al.ll2006h each are consistent with no vari- 
ation in c and with the gentle decline with increasing M 
expected for CDM models. In each case the scatter in 
log 10 c about the mean relation is not well constrained 
and is consistent with the CDM prediction. 

Our present investigation aims to improve significantly 
the constraints on the c — M relation by analyzing a 
wider mass range with many more systems than the 
previous Chandra and XMM studies. For this pur- 
pose we undertook a program to obtain accurate mass 
constraints on relaxed systems with 10 12 < M < 
10 14 M Q , representing 24 individual early-type galaxies 
and galaxy groups/clusters, each of whi ch possesses high- 
quali t y Chandra and/o r XMM data (|Humphrev et al.l 
120061 : iGastaldello et al.l 120061 : IZappacosta et all 120061 ) 
To these systems we added results for 1 5 rel axed, 
massive cluster s from iPointecouteau et ail (|2005l ) and 
IVikhlinin et al.l (|2006h . Using this combined sample, op- 
timized for X-ray mass measurements, we obtain empir- 
ical constraints on the local c — M relation with a simple 
power-law model and compare this model to the predic- 
tions of a suite of CDM models. 

We note that the primary purpose of this paper is to 
assemble accurate c — M measurements from Chandra 
and XMM observations from galaxy to cluster scales and 
to test whether c decreases with increasing M similar to 
that predicted by CDM models. However, we also make 
an initial attempt to constrain the parameters of CDM 
models to examine the efficacy of the c — M relation in 
this regard and to highlight sources of systematic error 
both in the models and data that need to be addressed 
with future studies. 

The paper is organized as follows. We describe the 
X-ray c and M measurements in $5] Our procedure for 
fitting a power-law model to the c—M data is discussed in 
IJU We define the cosmological models that are compared 
to the X-ray data in $H We present our results in $5] and 
discuss sources of systematic error in fj6] In fj7] we give 
our conclusions. All distance-related quantities, unless 
stated otherwise, are computed assuming a flat universe 
with Q m = 0.3, n A = 0.7, and H = lOO/i km s" 1 Mpc" 1 
with h = 0.7. 

2. OBSERVATIONS 



The measurements of concentrations and virial 
masses of 23 ea rly-type galaxies and galax y groups 
(|Humphrev et alJl2006tTGastaldello et al.ll2006l) form the 
core of our sample. In these papers we considered only 
those systems with the highest quality Chandra and/or 
XMM data. Moreover, to insure that hydrostatic equi- 
librium is a good approximation, we visually inspected 
all the early-type galaxies and groups in the public data 
archives and selected those systems that possessed the 
most regularly shaped X-ray images devoid of strong 
asymmetries. 

To these we added several clusters to populate the 
high-mass portion of the c — M diagram. Firstly, we 
include our analysis of t he radio-quiet (and ver y sym- 
metrical) cluster A2589 (|Zappacosta et ail l2006f l . Sec- 
ondly, we added c l usters from the XMM study by 
I Pointecouteau et "all (20051 ) and the Chandra study by 
Vikhl inTn et alj ( 2006 ) , each of whom also focused on 
the highest quality observations of the most relaxed sys- 
tems. For clusters t hat are common to b oth studies, 
we used results from IVikhlinin et al.l (|2006h because of 
more accurate background subtraction and temperature 
constraints, the latter resulting from the much smaller 
point spread function of Chandra; i.e., less biased tem- 
perature measurements in the presence of strong tem- 
peratu re gradients. For those few (lowes t mass) systems 
in the IPointecouteau et ail (|2005l ) and [Vikhlini n et al.l 
(2006) studies that overlap with our work cited above, 
we use our values. 

The final sample (see Table [lj consists of 39 systems 
spanning the mass range (0.06 — 20) x 10 14 M Q , with c 
and M values inferred using the NFW dark matter pro- 
file. We follow the convention that defines the virial ra- 
dius (r v i r ) so that the mean density within r v ; r equals 
Ap c , where p c (z) is the critical density of the universe 
at redshift z, and A is obtained from the solution to the 
top-hat spherical collapse mod el. To evaluate A we us e 
the approximation obtained bv lBrvan fc Norman! (|1998|) . 
For simplicity, when considering dark energy (DE) mod- 
els with w 7^ —1, we adopt the value of A appropriate 
for w = — 1. Although technically the viria l overdensity 
changes with w =/= — 1 (Ku hlen et al.|[2005T ). this choice 
has no effect on our conclusions. The values of A for 
the concordance cosmology at the appropriate redshift 
for each system are listed in Table [TJ 

For results obtained using A values different from those 
in Table [lj we co nverted c and M to ou r adopted A us- 
ing the formula of lHu fc Kravtsovl (|2003l ). which assumes 
an NFW mass profile. We converted the standard devia- 
tions for c and M in one of two ways. For those systems 
taken from our previous work, we have at least 20 Monte 
Carlo error simulations. Hence, we convert each simu- 
lated value of c and M to the new A and then compute 
the standard deviation of the converted values. This pro- 
cedure provides a self-consistent conversion of the stan- 
dard deviations. For systems taken from the literature, 
we simply convert the lower limit (e.g., c — er c ) and up- 
per limit (e.g., c + a c ) each to the new A. The new 
standard deviation is then set to one-half the difference 
of the converted upper and lower limits. (By perform- 
ing this procedure on our own data and comparing to 
the rigorous method using the error simulations, we find 
that this procedure is adequate for our present study.) 

Since we perform our analysis of the c — M relation 
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TABLE 1 

Halo Concentration and Virial Mass 
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Clog en 


(M e ) 


(M ) 


flog M 


Corr 


Corriog 


Ref. 
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(2) 


(3) 


(4) 


(5) 






(7) 


(8) 






(10) 


(11) 


(12) 


NGC4125 





0045 


101.5 


10.19 


1 


78 


7 


52c-02 


6.24c+12 


8.04c+ll 


6 


24c-02 


-7.87c-01 


-7.61c-01 




NGC720 


o 


0058 


101.6 


18.54 


7 


25 


1 


65e-01 


6.65c+12 


1.42c+12 


9 


63c-02 


-7.88c-01 


-7.83c-01 




NGC6482 


o 


0131 


102.3 


18.40 


5 


39 


1 


35e-01 


7.15e+12 


1.72e+12 


8 


63c-02 


-8.46c-01 


-8.87c-01 




NGC5129 


() 


0230 


103.1 


14.95 


2 


46 


7 


71e-02 


1.54e+13 


2.18e+12 


5 


96c-02 


-8.99g-01 


-8.94c-01 


b 


NGC1407 





0059 


101.7 


17.88 


4 


93 


1 


09c-01 


1.59e+l3 


3.12e+12 


9 


35c-02 


-8.41e-01 


-8.35e-01 




NGC533 


() 


0185 


102.7 


17.32 


1 


21 


3 


14c-02 


2.29e+13 


l.lle+12 


2 


05c-02 


-9.03c-01 


-9.12c-01 


b 


NGC2563 





0149 


102.4 


10.16 


3 


43 


1 


16c-01 


2.61e+13 


5.46e+12 


9 


78c-02 


-9.01e-01 


-9.74e-01 


b 


RGH80 


() 


0379 


104.4 


10.28 


1 


01 


4 


15c-02 


2.81e+13 


1.60c+12 


2 


55c-02 


-9.23c-01 


-9.28c-01 


b 


NGC1550 


() 


0124 


102.2 


17.16 


1 


12 


2 


80c-02 


3 12e+13 


1.49c+12 


2 


08c-02 


-9.53o-01 


-9.58e-01 


b 


NGC4472 


o 


0033 


101.4 


13.07 


1 


64 


5 


09c-02 


3 27p-t-1 3 




5 


97g-02 


-8 90p-01 


-8 89p-01 




NGC4325 


o 


0257 


103.3 


11.46 


1 


28 


5 


61c-02 




8.10e+12 


8 


34c-02 


-9.71c-01 


-9 83p-01 


b 


NOC4649 


o 


0037 


101.5 


20.80 


2 


45 


5 


17c-02 


3.50e+13 


4.62c+12 


6 


29c-02 


-8.88c-01 


-8.63c-01 




NGC5044 


o 


0090 


101.9 


11.32 





29 


1 


09c-02 


3.73e+13 


1.17e+12 


1 


40c-02 


-9.50c-01 


-9.53c-01 


b 


IC1860 


o 


0223 


103.1 


9.59 





88 


3 


91c-02 


5.10e+13 


4.62c+12 


3 


98c-02 


-9.66c-01 


-9.72c-01 


b 


MKW4 


o 


0200 


102.9 


12.72 





86 


2 


96c-02 


6.22c+13 


3.38e+12 


2 


39c-02 


-9.43c-01 


-9.46c-01 


b 


NGC4261 


o 


0075 


101.8 


3.72 





86 


1 


09c-01 


6.69e+13 


1.03e+13 


6 


35c-02 


-8.29c-01 


-7.95c-01 




RXJ1159. 8+5531 





0810 


108.0 


11.48 


2 


86 


1 


28c-01 


9.10e+13 


6.89e+13 


1 


73c-01 


-8.12e-01 


-9.79c-01 


b 


A262 


o 


0163 


102.5 


8.87 





67 


3 


46c-02 


l.lle+14 


1.06e+13 


3 


99c-02 


-9.77c-01 


-9.83c-01 


b 


MS0116. 3-0115 





0452 


105.0 


6.57 


2 


21 


1 


41e-01 


1.28e+14 


8.04e+13 


2 


14e-01 


-8.35c-01 


-9.72e-01 


b 


ESO5520200 


o 


0314 


103.8 


7.79 





75 


4 


46c-02 


1.31e+14 


1.96e+13 


6 


06c-02 


-9.53c-01 


-9.57e-01 


b 


MKW9 


() 


0498 


105.4 


7.48 





93 


5 


41c-02 


1.44e+14 


3.61e+13 


1 


llc-01 


-6.49c-01 


-9.36c-01 


d 


AWM4 


() 


0317 


103.9 


9.27 





84 


3 


72c-02 


1.62e+14 


1.82c+13 


5 


10c-02 


-9.18c-01 


-9.33c-01 


b 


ESO3060170 





0358 


104.2 


9.17 


1 


04 


fi 


12c-02 


1.81e+14 


7.56e+13 


1 


09c-01 


-8.84e-01 


-9.76e-01 


b 


A2717 


o 


0490 


105.3 


6.45 





36 


2 


46c-02 


1.83e+14 


1.20c+13 


2 


84c-02 


-9.79c-01 


-9.83c-01 


b 


A1991 


o 


0592 


106.2 


8.91 


o 


66 


3 


22c-02 


1.93e+14 


2.67e+13 


6 


04c-02 


-l.OOc+00 


-1.00e+00 




A1983 





0442 


104.9 


5.34 





09 


g 


15c-02 


1.99c+14 


7.62c+13 


i 


76c-01 


-2.87c-01 


-3.92c-01 


d 


A2589 


o 


0414 


104.7 


6.36 





32 


2 


33e-02 


3 24e+14 


2 37o-)-1 3 


3 


00e-02 


-9 5fip-01 

U • UUt V;l 


-9 6fip-01 




A2597 





0852 


108.4 


8.26 





70 


3 


71e-02 


3.56e+14 


3.91c+13 


4 


80e-02 


-7.87e-01 


-l.OOc+00 


d 


A383 





1883 


116.8 


9.61 





89 


4 


03e-02 


4.70c+14 


4.76c+13 


4 


41e-02 


-5.12e-01 


-l.OOc+00 


e 


A133 





0569 


106.0 


6.72 





61 


3 


97e-02 


5.37c+14 


6.44c+13 


5 


23e-02 


-5.50e-01 


-l.OOc+00 


c 


A1068 





1375 


112.7 


5.44 





38 


3 


07e-02 


6.96e+14 


6.01c+13 


3 


76e-02 


-9.42e-01 


-l.OOc+00 


d 


A907 





1603 


114.5 


7.75 





94 


5 


27e-02 


7.39e+14 


5.99e+13 


3 


53e-02 


-3.87e-01 


-l.OOc+00 


e 


A1795 





0622 


106.5 


6.80 





38 


2 


44e-02 


1.02e+15 


8.78c+13 


3 


75e-02 


-6.48e-01 


-l.OOc+00 


c 


PKS0745 





1028 


109.9 


7.32 





57 


3 


40e-02 


1.20e+15 


1.44c+14 


5 


24e-02 


-2.64e-01 


-l.OOc+00 


d 


A478 





0881 


108.6 


7.61 





58 


3 


29e-02 


1.25e+15 


1.65e+14 


5 


74e-02 


-2.28e-01 


-l.OOc+00 


e 


A2029 





0779 


107.8 


8.47 





44 


2 


26e-02 


1.27c+15 


1.18c+14 


4 


02e-02 


-4.19e-01 


-l.OOc+00 


e 


A1413 





1429 


113.1 


6.56 





38 


2 


52e-02 


1.29e+15 


1.29e+14 


4 


37e-02 


-4.41c-01 


-l.OOc+00 


c 


A2204 





1523 


113.9 


6.77 





55 


3 


51c-02 


1.41c+15 


1.55e+14 


4 


80e-02 


-2.55e-01 


-l.OOc+00 


d 


A2390 





2302 


120.0 


4.13 





32 


3 


41c-02 


2.13c+15 


2.14c+14 


4 


38o-02 


-3.13o-01 


-l.OOc+00 


e 



Note. — Col.(l): Cluster name. Col. (2) Rcdshift. Col. (3) Reference overdensity for virial radius and mass definition. Col. (4) Scaled concentration 
Co — (1 + z)c. Col. (5) Standard deviation on cq. Col. (6) Standard deviation on log 10 Co. Col. (7) Virial mass. Col. (8) Standard deviation on M. 
Col. (9) Standard deviation on log 10 M. Col. (10) Correlation co effic ient for cq and M . Col .(ll) Correlation coefficient for l og 1n CQ and log 1n M. 
Col.(1 2) R eference for c and M va lues: a — Humphrey ct al. ( 2006), b - Gastaldcllo ct al. ( 20 061) , c — IZap pacosta ct al. ( 2006) , d — Pointccoutca u et al.l 
(2005), c - Vikhlinin ct al. (2006). Note that literature results obtained at different A values have been converted to those in column (3) as explained 
in 



in log space (see below), we also require the standard 
deviations of log 10 c and log 10 M. As above, for sys- 
tems analyzed in our previous studies, we compute self- 
consistently the log errors using Monte Carlo error sim- 
ulations. The other systems use the simple prescription 
noted above, but using the log values. The log errors on 
c and M are listed in Table Q] as, respectively, criogc and 

OlogM- 

For the systems with Monte Carlo error simulations 
we also computed the covariance between c and M 
(and log 10 c and log 10 M) for the different systems; e.g., 
XX C * — ( c ))(^i ~ (M)), where the brackets represent the 
mean quantity, and the sum is over the simulations. In 
Table Q] we report these as the correlation coefficient, 
which is the covariance divided by a c au (or 0i ogc ciogM)- 
For the objects from the literature we used the average 
covariance obtained from the others to evaluate the cor- 
relation coefficient, setting the coefficient to -1 if a value 
< — 1 was obtained. Most of the correlation coefficients 
have values near -1 indicating that c and M are anti- 
correlated. 

Finally, the results listed in Table Q] for the concen- 



tration parameter actually refer to the quantity, cq = 
(l + z)c, since this quantity is what is fitted in our analy- 
sis below. Note that Co is also used when computing the 
covariance and correlation coefficient. 

3. ANALYSIS METHOD 

We focus our analysis on a simple power-law represen- 
tation of the c — M data, 



where z is the redshift, and both C14 and a are con- 
stants independent of M . We set the reference mass 
to, M14 = 10 14 /i _1 M©, which lies close to the midpoint 
(in log space) of the mass range of our sample. CDM 
models generally predi ct that c decreases with increas- 
ing M (e.g., B01; D04: iKuhlen et al.ll2005t iMaccio et all 
2006§. A key goal of our study, therefore, is to deter- 
mine whether a < 0. For the most massive systems 
(M > 3 x 1O 14 M ) in CDM models D04 hnd that c u 
is much more sensitive to cosmological parameter varia- 
tions than a (see Q. 
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TABLE 2 

COSMOLOGICAL MODEL PARAMETERS 



Name 






U B h 2 


h 


CT8 


n s 


-w 


ACDM1 


0.30 


0.70 


0.022 


0.7 


0.90 


1.00 


1.0 


ACDM3 


0.24 


0.76 


0.022 


0.73 


0.76 


0.96 


1.0 


DECDM 


0.30 


0.70 


0.022 


0.7 


0.90 


1.00 


0.6 


QCDM 


0.30 


0.70 


0.022 


0.7 


0.82 


1.00 


as 0.8 


OCDM 


0.30 


0.00 


0.022 


0.7 


0.90 


1.00 





Note. — f2 m is the energy density parameter for matter in the 
universe; Q\ is the energy density parameter associated with a 
cosmological constant or, more generally, dark energy; Clg is the 
energy density parameter of baryons; h is f/o/100 km s _1 Mpc - 1 ; 
erg is the rms mass fluctuation within spheres of comoving radius 
8k' 1 Mpc. See gU 

To constrain the parameters of eqn. [Tj it is neces- 
sary to account for the error estimates on both c and 
M. Consequently, we employ the BCES method (i.e., bi- 
var iate correlated errors with i ntrinsic scatter) described 
by lAkritas fe Bershadvl (|1996f ) to estimate a, cr 4 , and 
the intrinsic scatter about the best relation. We per- 
formed the BCES fitting using software kindly provided 
by M. Bershady 4 . Since BCES is a linear regression 
method we transform eqn. [JJ to the form, y = ax + b, 
where x = \og 1Q M, y = log 10 c = log 10 (l + z)c, and 
b = log 10 (cu/Mf^). The parameter en is derived from 
a and b: cu = 10 b M{* 4 . Since the fractional error 
flog m / log 10 M is typically smaller by more than an or- 
der of magnitude of the fractional error, <7i og Co / log 10 Co 
(see Table [JJ), we always use x = log 10 M as the indepen- 
dent variable. By default we also use the Corri og values 
in Table [TJ for the BCES method. Despite the strong 
anti-correlation between c and M, we find that includ- 
ing the Corriog values in the analysis has an insignificant 
impact on the estimated a and cu values. 

We determine the best estimate of a and C14 by per- 
forming the BCES method on 10 6 bootstrap resamplings 
of the data. We take the mean of the bootstrap simu- 
lations to be the best estimates and construct error el- 
lipses about these best values. The confidence contour 
spacings computed from the bootstrap simulations corre- 
spond quite closely to the A^ 2 values for two parameters 
assuming normally distributed errors. 

We estimate the intrinsic scatter on the concentration 
as, 

(°y)int ~~ ( CT a) total _ (°Wstat 

N 1 N 

^("^-i^- (2) 

i=l i=l 

where N is the number of data points (xi,yi) = 
(log 10 Mi, log 10 co.i) corresponding to the entries in Table 
[TJ y™- odel = axi+b, where a and b are the best (mean) es- 
timates from the bootstrap simulations; and a Vi = <J\ ogCo 
in Table [[J 

4. COSMOLOGICAL MODELS 

In Tablc[2]wc define the CDM-bascd cosmological mod- 
els to be compared with the observations. ACDM1 is the 
standard concordance model, which essentially reflects 
the combined constraints from the first year of WMAP 

4 http:/ /www. astro. wisc.edu/~mab/archive/stats/stats. html 



CMB observations, the supernovae Hubble diagram, and 
the l arge-scale clustering of galaxies (e.g.. ISpergel et ail 
120031 ). ACDM3 effectiv ely updates ACDM 1 using the 
third year of WMAP (jSpergel et all l2006f) . Each of 
these models assumes a constant dark energy equa- 
tion of state with w — —1.0 (i.e., a cosmological con- 
stant). In addition, we consider a dark energy model with 
w = -0.6 (DECDM) and a quint essence model (QCDM) 
with a Ratra-Peebles potential (jRatra fe Peebles! 119881 : 
iPeebles fc Ratrall2003l and references therein) as imple- 
mented in D04. The QCDM model has a nearly constant 
w(z) pa —0.8 (see Fig. 1 of D04) and a lower er s . Finally, 
we include an open model (OCDM) with the same pa- 
rameters as ACDM1 except with fl\ = 0. 

The median relation between concentration and virial 
mass as a function of redshift for CDM halos is described 
well by the semi-analytic model proposed by B01, 

c(M,z) = k\±^, (3) 

14-2 

where K is the normalization constant, z is the redshift 
of the halo, and z c is the redshift when the halo "col- 
lapsed" . This collapse redshift is defined implicitly by 
the equation, 

a(FM) = 5 c (z c ), (4) 

where 6 C is the equivalent linear overdensity for spherical 
collapse at z c , a is the z — linear rms density fluctua- 
tion corresponding to a mass FM, and F is a constant. 
(Note that eqn. [3] is not a power-law, though over small 
mass ranges - within a factor of 5-10 - the B01 model 
can be well approximated by one.) The constants F and 
K must be specified by comparison with numerical sim- 
ulations. 

For halos with M < 10 13 /i _1 Af Q it is found that 
F = 0.01 and K w 3.5 provide the best description of 
the results from N-body simulations, although F = 0.001 
and K » 3 is a cceptable (B01; iKuhlen et all 120051 : 
iMaccio et al.l 12006V Larger box simulations that focus 
on higher mass ranges clearly prefer F = 0.001. One of us 
(J. Bullock) finds that F = 0.001 and K = 3.1 matches 
the c-M res ults for CDM halos (0.6 -2.5 x IO^/^Mq) 
simulated bvrTasitsiom i et all (|2004T ). For simulations of 
even more massive clusters (3.1 — 17 x 10 14 /i~ 1 M Q ), D04 
find higher concentrations and that K — 3.5 is required 
for the same F. Consequently, when comparing CDM 
models to the full mass range of the X-ray data in Table 
[1] we shall consider the c — M relations obtained with 
the B01 model using F = 0.001 with both K = 3.1 and 
K = 3.5. (In £)6.3I we discuss further the different F 
and K values obtained by different investigators.) In the 
high-mass "cluster" regime D04 find that the c — M re- 
lation for CDM models is adequately parameterized by 
the simple power-law model given by eqn. [TJ. We list 
the power-law parameters obtained by D04 in Table [3] 
but converted to our definition of virial radius. (D04 
use A = 200p b = 60p c for Q m = 0.3.) That is, ini- 
tially we compute c(M) using C14 and a from D04 (at 
z = 0). Then we convert the resulting c and M values 
to A = 101.1 using the approximation of lHu fc Kravtsovl 
(|2003h . These values are used to compute the slope and 
C14, which are listed in Table [3l The parameter C14 de- 
creases by w 20% while a remains nearly constant under 
this transformation. For comparison, in Table [3] we also 
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TABLE 3 

Power-Law Approximation for CDM Clusters 



Name 


a 


C14 


ACDM1 


-0.104 


7.63 


DECDM 


-0.094 


9.05 


QCDM 


-0.111 


7.39 


OCDM 


-0.091 


11.47 


ACDM1 B 


-0.142 


8.43 


ACDM3 b 


-0.155 


6.23 



5. RESULTS 

In [JB] we discuss possible sources of systematic error 
both in the measurements from the X-ray data and the 
theoretical predictions. Possible systematic errors affect- 
ing the X-ray mass measurements of our sample of re- 
laxed galaxies and clusters are expected to be small and 
do not change the conclusions of the analysis presented 
here (see beginning of fj6] and j j6.2| . 



Note. — Power-law parameters (eqn. [TJ for the theoretical re- 
lation between concentration and virial mass of high-mass clusters 
obtained by D04 for the ACDM1, DECDM, QCDM, and OCDM 
models. We converted the results presented in Table 2 of D04 to 
our definition of the virial radius at z = (i.e., A = 101.1). The 
ACDMIb and ACDM3b models are the same as the ACDM1 and 
ACDM3 models, but we computed the power-law parameters using 
the B01 model with F = 0.001 and K = 3.5. See 



list a and cu computed from the tangent line to the 
c - M profiles located at M = M u of the ACDM1 and 
ACDM3 models obtained using the B01 approach with 
F = 0.001 and K = 3.5. These results are listed as 
ACDMIb and ACDM3 B respectively. (We note that the 
virial quantities for the X-ray data in Table [1] all refer to 
the ACDM1 cosmology. For consistent comparison, we 
also use the ACDM1 cosmology when converting all the 
D04 models mentioned above.) 

The concentration is sensitive to the f ormation time 
and dynamical state of the halo (e.g.. iNavarro et~al 



1997; Eke et all 120011: Uind 120001 : B01; IWechsler et al 
20021: IMaccio et al 1 120061) . For example, Uind (|2000f ) 
and Maccio et al.l (|2006( ) find that their "relaxed" dark 



halo samples have concentrations that are systematically 
larger by ~ 10% co mpared to the whole population. 
IWechsler et al.l ()2002l ) found that the halo concentration 
at fixed mass is set almost exclusively by the "formation 
epoch" of the halo (the time when the mass accretion 
rate of the halos is slowed below a critical value). In- 
deed, when IWechsler et~ail (|2002j) focused on halos with 
no major mergers since z = 2 they found ~ 10% higher 
concentrations for that population. Therefore, a proper 
comparison between the theoretical c — M relation with 
observations requires the observed and simulated halos 
be selected in a consistent manner. The X-ray data pre- 
sented in Table Q] represent objects that were selected 
to be the most relaxed, X-ray-bright, systems. But the 
theoretical c — M models discussed above (Table [3]) were 
obtained using all available halos in the N-body simula- 
tions. Consequently, we consider this source of system- 
atic error in our comparisons of the CDM models with 
the X-ray data. 

Finally, the magnitude of the intrinsic scatter about 
the median theoretical c—M relation does not vary over a 
large class of CDM models and is independent of M. The 
value of oiogc « 0.14 (log = log 10 ) obtained by B0 1 has 
been found by several indep endent investigations dJind 
120001: iTasitsiomi et~aLl 120041: D04; IMaccio et all [2006); 
note this consistency is observed after accounting for the 
fact that In, rather than log 10 , is used in these other 
studies. Early forming, rela xed, halos tend to exhibit 
smaller scatter ai g C ~ .10 (pind l200"ol ; IWechsler et al.1 
2002: IMaccio et al.ll20"0fil) . 



5.1. All Halos 

When fitting the power-law relation (eqn. [T]) to the 
entire sample we obtain a = —0.172 ± 0.026 and C14 = 
9.0 ± 0.4 (quoted errors la). The estimated intrinsic 
scatter in log 10 (1 + z)c is 0.102 ± 0.004 (eqn. |2J), with 
a total scatter of 0.12. In Figure Q] we plot the best- 
fitting model with intrinsic scatter. Inspection of Figure 
[T] reveals that the power-law with intrinsic scatter is a 
good approximate representation of the X-ray results, 
especially for log 10 M > 13.5. For lower masses the 
X-ray measurements tend to lie above the best-fitting 
model, though most are consistent with lying within the 
lcr ra nge covered by the intri nsic scatter. Using equation 
[3] of iTremaine et al.l (|2002f l we compute a reduced % 2 
value of 1.07 (37 dof) for the fit, confirming the visual 
impression that the power law is a good, but not perfect, 
representation of the X-ray data. 

The slope a is constrained to be negative at the 6.6c 
level, demonstrating at high significance that c decreases 
with increasing M, as expected in CDM models. In Fig- 
ure Q] we also plot the theoretical c — M relation for 
the ACDM1 model for both K = 3.1 and K = 3.5 
as discussed in $4j For log 10 M > 14 the K = 3.5 
model recommended by D04 for clusters lies very close 
to, but just below, the empirical power-law fit. Allow- 
ing for a systematic rs 10% in crease in concentration for 
relaxed, early forming h alos (jJingl 120001 : IWechsler et al.1 
l2002UMacci6 et al.ll2006l ) lifts the K = 3.5 model slightly 
over the power-law, with a marginal improvement in 
the agree ment between the models. A recent theoreti- 
cal study (jWang et al.l 12006) finds that the higher mass 
halos experience less of bias associated with relaxation 
and formation time. In fact, systematically increasing 
the concentrations of the K = 3.5 model by ~ 7% would 
provide a closer match to the empirical power-law fit to 
the X-ray data. 

For log 10 M < 14 the K = 3.5 model must quickly 
transform into the K = 3.1 model as M decreases from 
the cluster to the galaxy/group regime as discussed in 
21 But in this mass range the K = 3.1 model lies con- 
sistently below the power-law model and is similar to 
the la lower limit on the power-law given by the intrin- 
sic scatter. Increasing the concentrations of the K = 3.1 
model by w 10% for relaxed halos improves the agree- 
ment with the power-law, and has the effect of essentially 
pushing the K = 3.1 model into the K = 3.5 model. 
With this accounting for the systematic bias arising from 
relaxed, early forming halos the ACDM1 model is an 
equally good representation of the X-ray data for masses 
log 10 M > 13.5. For the lowest masses, the ACDM1 
model still lies below the power-law (within the la in- 
trinsic scatter) and may represent a real discrepancy. 

We summarize this qualitative discussion with a quan- 
titative comparison of the ACDM1 model with the data 
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Fig. 1. — Results of fitting eqn.[TJin log space to the entire sample. Displayed are the best-fitting model (solid blue line), the la intrinsic 
scatter (dotted blue lines), and the predicted relation obtained from cosmological simulations for the (Left panel) ACDM1 and (Right panel) 
ACDM3 models (dashed lines). The two dashed lines attempt to represent fits to different mass ranges in the cosmological simulations as 
explained in the text (Q. The lower (red) dashed line refers to fits of halos up to M ~ 0.3 X fO 15 M0 (K = 3.f) while the upper (green) 
dashed line refers to fits of higher mass halos M » (0.3 - 1) X f0 15 M Q (K = 3.5) obtained by D04. 
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Fig. 2. — Best-fitting value (X) and confidence contours (68%, 95%, 99%, 99.9%, 99.99%) for a and C14 obtained from the entire sample. 
(Left paneJ) As in Fig. [TJ results for the ACDMf (green i) and ACDM3 (red 3) models are displayed. Models with subscript "c" refer to 
the higher normalization K = 3.5 for t he cluster regime pr eferred by D04. (Right panel) Models have C14 increased by i0% to represent 
relaxed, early- forming halos (see [|4]and Maccio ct al. 2OQ3). 



as represented by the results of the empirical power-law 
fit. (We do not fit the B01 model to the data for the 
following reasons. First, it is non-trivial to account for 
the error bars in both coordinates when fitting a gen- 
eral model; the Akritas & Bershady method only ap- 
plies to linear regression. Second, the prediction of the 
B01 model possesses significant systematic uncertainty 
in the parameter K (at least for masses < 10 14 M Q ) as 
discussed above. By showing the B01 model predictions 



separately for the range of interesting K values, without 
fitting, we provide a clear demonstration of the impor- 
tance of reducing this systematic uncertainty with future 
theoretical studies. When that happens, formal fitting 
of the B01 model will be investigated.) In Figure [2] we 
show the error contours estimated for a and C14 from 
the bootstrap simulations (SJ3J). We calculated the slope 
and concentration of the c — M relation of the ACDM1 
model at lO^h^Mp, for both the K = 3.1 and K = 3.5 
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cases. (The values for K = 3.5 are listed in Table [3l) 
These results are plotted in Figure [2] as "1" for K = 3.1 
and l c for K = 3.5, the latter representing the "cluster" 
regime. The right panel of Figure [2] increases the con- 
centrations by 10% to represent relaxed, early forming 
systems, and shows the ACDM1 model lies near the 68% 
confidence contour. This level of agreement applies for 
masses log 10 M > 13.5. As discussed above, the agree- 
ment is worse at the lowest masses. However, perform- 
ing a similar comparison at 10 13 /i _1 M we find the local 
ACDM1 slope and normalization lie just within the 95% 
confidence contour of the power-law model (both with 
and without the 10% correction for relaxed systems); 
i.e., the disagreement is not very significant even at the 
low-mass end. If real, the discrepancy may signify very 
early forming fossil groups in the sample (D'Onghi a et al.l 
f2005h . 

The intrinsic scatter of 0.10 in log 10 (1 + z)c obtained 
for the power-law fit is smaller than the value of 0.14 
for all dark matter halos in CDM simulations but agrees 
extremely well with the value expected f or the most re- 
laxed , early forming syst ems (|Jing||2000t IWechsler et al.1 
l2002t iMaccio et al.1 l2006h . Consideration of this result 
for the scatter and the average c — M relation above, we 
conclude that the ACDM1 model is consistent with the 
X-ray data, provided the X-ray sample reflects the most 
relaxed, early forming systems in the population. (This 
corroborates our selection criteria discussed in ij^l) 

Now we perform the analogous comparison of the 
ACDM3 model with the X-ray data and associated 
power-law fit. The c — M relation is displayed in the 
right panel of Figure [TJ The ACDM3 model lies well- 
below the power-law at all masses, even when allowing 
for the expected 10% increase in concentration for re- 
laxed halos. In Figure we compare the ACDM3 model 
at 10 1 K~ x Mq with the empirical power-law. Even when 
considering the 10% increase in concentrations for re- 
laxed halos, the ACDM3 model lies on the 99.99% con- 
tour. Systematic errors associated with the X-ray mea- 
surements cannot explain this level of disagreement (see 
ffl- 

The key parameter responsible for the poor perfor- 
mance of the ACDM3 model with respect to ACDM1 is 
the low value of as (0.76). In order to bring the ACDM3 
model within the 99% contour in the right panel of Fig- 
ure [2] requires as > 0.84 where we have kept the other 
cosmological parameters fixed to their values in Table 
[21 This limit is conservative since (1) we use K = 3.5, 
(2) we assume a full 10% upward shift for the bias from 
relaxed, early forming systems, which may be less for 
massive clusters, and (3) we have approximated the B01 
models as power-laws using their predictions only near 
1O 14 /i _1 M0. The sensitivity of the concentrations to as 
results from the impa ct that as has on the average halo 
formation times (e.g.. [Eke et al.ll200ll lAlam et al.l 120021 : 
Ivan den Bosch et al.ll2003[) . 

Other cosmological parameters, however, contribute to 
the large discrepancy of the ACDM3 model. If we re- 
move the tilt of the power spectrum (i.e., set n s = 1) 
then we obtain as > 0.80 at the 99% confidence level. 
Finally, if we further set f2 m = 0.3 and h = 0.73 (so 
the model is the same as ACDM1 except with variable 
as), then the 99% constraint falls to as > 0.76, the value 
associated with the ACDM3 model. (We reiterate this 




log,o M 

Fig. 3. — Results of fitting eqn.[T]m log space only for halos with 
M > 10 14 Mq. Displayed are the best-fitting model (solid blue 
line) and the ltr intrinsic scatter (dotted blue lines). Also shown 
is the predicted relation for the ACDM1 model from cosmological 
simulations. The (solid red) line is the power-law fit obtained by 
D04 while the (dashed red) line is the B01 model with F = 0.001 
and K = 3.5 found by D04 to best match the simulated clusters. 
Note the D04 model is converted to our definition of the virial 
radius (Q. 

limit is conservative as noted above. Similarly, we obtain 
as < 1.07 at 99% confidence, where the K = 3.1 model is 
used here to be conservative.) Hence, although the lower 
value of as is the primary cause of the poor performance 
of the ACDM3 model, the combined action of the power- 
spectrum tilt with the lower value of f2 m exacerbate the 
discrepancy with the X-ray c — M relation. 

If the ACDM3 parameters are correct, particularly the 
low value of as, then a fundamental modification of the 
model is required to increase the concentra tion values to 
match the X-ray results. Both D04 and iKuhlen et al.l 
(2005) have shown that changing the dark energy equa- 
tion of state parameter (w) has the effect of systemati- 
cally raising (larger w) or lowering (smaller w) halo con- 
centrations. As we show below, a model with w w —0.8 
and as ~ 0.8 can describe the X-ray data in the cluster 
regime. 

5.2. Only Halos with M > 10 14 M Q 

As discussed in SJU D04 provide results of fitting eqn. 
[I] to a suite of CDM models for massive clusters, in- 
cluding an open model and multiple dark energy models. 
We may compare these theoretical predictions directly 
to the power-law fit to the X-ray data with the following 
considerations. Firstly, we convert D04's results to our 
definition of the virial radius (see S}4]). Secondly, since 
D04 analyze only massive clusters (M > 4 x 10 14 M Q , 
h = 0.7 and converted to A = 101.1), we restrict our 
analysis to high-mass systems as well. In order to allow 
more precise constraints on the power-law fit we con- 
sider clusters down to a somewhat smaller mass limit, 
Al > 1 x 10 14 Mq. (As shown below, this choice is jus- 
tified since there is no obvious trend in the c — M rela- 





M4 

Fig. 4. — (Left panel) c-M relations for various CDM models obtained by D04 for galaxy clusters (see Table [3}. The solid (red) line is 
ACDM1, the dashed (blue) line is OCDM, the dotted (magenta) line is DECDM, and the dot-dashed (cyan) line is QCDM. Finally, the 
short dashed (green) line is the ACDM3 model computed using the B01 prescription with K = 3.5. (Right panel) Best-fitting value (X) 
and confidence contours (68%, 95%, 99%, 99.9%, 99.99%) for a and c i4 obtained by fitting only systems with M > 1O 14 M . Also shown 
arc the results for the models plotted in the left panel: ACDM1 (1), OCDM (O), DECDM (D), QCDM (Q), and ACDM3 (3 B ). 

models are very similar and provide a reasonable descrip- 
tion of the X-ray data, though the B01 representation is 
slightly steeper. 

In the left panel of Figure [4] we plot the power-laws 
obtained by D04 for the CDM models as well as the B01 
representation of the ACDM3 model with K = 3.5. The 
values of a and cu for all of these models (Table [3]) are 
plotted in the right panel of Figure [4] along with the 
error contours derived from the power-law fit to the X- 
ray data. As done in the previous section, the a and cu 
values for the B01 ACDM3 model were obtained using 
the c — M slope near M = lO 14 /^ 1 M Q . 

We emphasize that the D04 parameters may be com- 
pared directly to the results of our power-law fit to the 
X-ray data, since the same model is fitted over nearly 
the same mass range. Since the B01 representation of 
the ACDM3 model is very nearly a power-law over the 
range being investigated, a similarly direct comparison is 
appropriate. 

The results for the ACDM1 and ACDM3 models ob- 
tained using only the high-mass subsample corroborate 
those obtained from the full sample. As noted above, 
the ACDM1 model is an acceptable match to the X-ray 
data; the parameters lie near the 68% confidence con- 
tour in Figure [4j even if c is increased to correct for 
early forming halos. In contrast, the ACDM3 model is 
very inconsistent with the high mass data, with parame- 
ters lying outside the 99.99% contour. Even if we apply 
the full 10% correction to Cu appropriate for the most 
relaxed, early forming systems, then the model merely 
moves on top of the 99.99% contour. However, as noted 
above, it is expected that this correction is less for the 
highest mass systems. Hence, the ACDM3 model is also 
rejected using only the most massive clusters. 

The c — M data for high-mass clusters clearly exclude 
the OCDM model. The concentrations lie systematically 



tion over this mass range.) Finally, we need to consider 
the biases (« 10% higher c, smaller scatter in log 10 c) 
for relaxed, early forming systems as done for the whole 
sample in the previous section. However, the systematic 
increase in c may be less since we are considering the 
most massive systems ([Wang et a l. 2006). 

In Figure [3] we show the result of fitting eqn. p] to the 
high-mass sample, for which we obtain the following pa- 
rameter constraints, a = — 0.04±0.05 and C14 = 7.6±0.5. 
These parameters are consistent with a constant c — M 
relation, as well as profiles that both decrease and in- 
crease with increasing M, and agree well with similar fits 
obtained in this m ass range from previous X-r ay studies 
with either XMM (jPointecouteau et al.ll2005[ ) or Chan- 
dra (jVikhlinin et alJ l2006f >. Since the power-law fit is 
quite consistent with a = 0, it follows that it is the lower 
mass systems (< 10 14 M Q ) which require a < found 
whe n fitting the whole sam ple in the previous section 
(see iGastaldello et all [20061) . (Note that the power-law 
parameters inferred in the previous section for the full 
sample lies on the 99% confidence contour obtained for 
the fit to the high-mass subsample - see Figure [H) 

We measure an intrinsic scatter in log 10 (1 + z)c of 
0.07 ±0.01 which is somewhat less than the value of 0.10 
obtained for the entire sample. D04 also obtain 0.10 for 
the intrinsic scatter (after accounting for a factor of 2.3 
converting In to log 10 ) f° r a ^ °f their CDM clusters. The 
smaller scatter we measure presumably can be explained 
because the X-ray sample preferentially contains relaxed, 
early forming systems. 

For comparison we plot two versions of the ACDM1 
model in Figure[3] The first version we show is the power- 
law fit obtained by D04 converted to our definition of the 
virial radius (see []4]). The second version arises from us- 
ing the B01 model with K = 3.5, which D04 found to 
best represent their high-mass cluster simulations. Both 
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above the data and the ACDM1 model. The latter is 
expected because structures form earlier in the OCDM 
model. In principle, the OCDM model may be brought 
into acceptable agreement with the X-ray data by lower- 
ing o~%. We estimate the effect of lowering <7g using the 
ACDM1 B01 model. For example, lowering erg from 0.90 
to 0.70 lowers c u by 25%. This would move the OCDM 
model to the right in Figure [4] on to the 99% confidence 
contour. Considering any of the 10% systematic bias 
expected for early forming systems would increase the 
discrepancy. (Note that for erg = 0.76, appropriate for 
the ACDM3 model, the OCDM model would be rejected 
at the w 99.9% confidence level without additional con- 
sideration of early formation bias.) 

Studies of the abundances of galaxy clusters using weak 
gr avit at i onal lensing typically find ag consistent with 
0.9 (e.g-.lHoekstra et al Jl2002t lVan Waerbeke et al.ll2005l; 
iJarvis et al.l l2003fh The optical study of iRines et all 
(|2006j l uses the caustic technique to measure cluster 
abundances and also finds <7g rs 0.9 with a lower limit 
of 0.72 (95% conf.). However, X-ray studies of clus- 
ter abundances show a large variation in erg , with val- 
ues ranging from 0.7-0.9 (see lArnaudl I2005L and refer- 
ences therein). The X-ray studies use a simple con- 
version between mass and global temperature or mass 
and global X-ray l uminosity (L x ). A recent study by 
iStanek et al.l ([2001 argues that considering reasonable 
errors in the M — L x conversion allows agreement of clus- 
ter abundances with the recent WMAP parameters wit h 
cr 8 = 0.85 and Q m = 0.24 - though see iReiprichl (|2006h . 
We take erg = 0.7 to be a conservative lower limit estab- 
lished by the abundances of galaxy clusters. 

Therefore, an open CDM model with f2 m rs 0.3 is ruled 
out (> 99% confidence) from joint consideration of only 
the X-ray c — M relation and cluster abundances. Since 
as and O m are tightly coupled from cluster abundance 
studies, and all the studies cited above give er 8 « 0.45 
for f2 m w 1, using the X-ray c — M data we find that we 
can also exclude CDM models with f2 m w 1 at high sig- 
nificance (> 99.9% confidence). For this comparison we 
have employed the entire sample of 39 systems and have 
set h w 0.5 for the fi m fss 1 models to better satisfy cos- 
mic age constraints from stellar populations in globular 
clusters (|Chaboyerlll998l ). Consequently, the combina- 
tion of the X-ray c — M relation and cluster abundances, 
assisted by constraints on the cosmic age, provides novel 
evidence for a flat, low-f2 m universe with dark energy 
using observations only in the local (z <C 1) universe. 

Finally, we consider the alternative dark energy mod- 
els, DECDM and QCDM. The DECDM model lies sys- 
tematically above the ACDM1 model and deviates from 
the power-law fit to the X-ray data at the 99.9% confi- 
dence level. This behavior is expected since w = —0.6 
implies halos form earli er compared to th e w = — 1 
ACDM1 model (see also iKuhlen et all 12005( 1 . By lower- 
ing as it should be possible to bring the DECDM model 
into good agreement with the data. To illustrate this 
trade-off between w and ag, we consider the QCDM 
model which has w Rs —0.8 and erg = 0.82. As shown 
in Figure [4] the QCDM model matches the X-ray results 
nearly as well as the ACDM1 model. In fact, a mild in- 
crease in C14 owing to effect from selecting relaxed, early 
forming halos would produce even better agreement. 

Our analysis indicates that a model similar to QCDM 



is able to satisfy the most recent constraints from other 
cosmological observations since w w —0.8 is marginally 
consistent with the constraints imposed b y supernovas 
(w = 1.023 ± 0.090 (stat) ± 0.054 (sys), lAstier et al.l 
2006), and erg = 0.82 is also marginally consistent with 
the 3-yr WMAP results. Put another way, if the low ag 
value of ACDM3 is correct, then the X-ray c — M data 
imply — 1 < w ~ —0.8. 

6. SYSTEMATIC ERRORS 

When comparing to theoretical predictions of the c—M 
relation we have considered the small biases in the mean 
concentration level and scatter expected for the highly 
relaxed, early forming systems selected for our study. We 
found that the ACDM1 model is a good representation 
of the X-ray data provided this bias for relaxed systems 
applies, especially for the low-mass range of our sample 
(M < 10 14 M©). However, the ACDM3 model does not 
agree with the c — M data, even considering this bias. 

Could other sources of systematic error seriously af- 
fect these conclusions? For 24 of 39 systems in Ta- 
ble Q] we have performed a d etailed investigat i on of 
systematic errors on c and M ([Humphrey et al.l 120061 : 
IZappacosta et al.l 120061 : iGastaldello et al J 12006( 1 Gener- 
ally, we find that the estimated systematic error is less 
than the statistical error, particularly for the systems of 
lowest mass. For the higher mass systems, the errors 
are usually comparable. But from consideration of all of 
these systems we recognize no obvious trend that would 
systematically shift the observed c — M relation in one 
direction. 

Although we do not believe systematic errors associ- 
ated with the X-ray data analysis compromise the con- 
clusions of our present investigation, here we list the most 
important issues to be resolved in future studies that 
seek to use the X-ray c — M relation for precision con- 
straints on cosmological models. Below in §6.21 we also 
consider the impact on our present analysis of systematic 
errors resulting from the hydrostatic equilibrium approx- 
imation. 

(Note we verified that our analysis is insensitive to the 
definition of the virial radius. An interesting example is 
to compare to results obtained for A = 2500, because 
within this radius all of the systems in our study possess 
good X-ray constraints. When defining the virial radius 
to correspond to A — 2500 for every system, we arrive 
at the same conclusions obtained in §21 but with corre- 
spondingly different a and cn values. The predictions of 
the cosmological models were also converted to A = 2500 
for this comparison.) 

6.1. Early Formation Bias 

Accounting for the biases associated with preferentially 
observing the most relaxed, early forming systems is re- 
quired for prec i sion constraints on cosmology. While 
IWechsler et all (|2002h found a nearly one-to-one corre- 
lation between halo "formation epoch" and concentra- 
tion, it is not obvious how this theoretically motivated 
parameter should connect with the dynamical state of 
a real cluster. More direct quantifiers (e.g. rejecting 
cases with recent major mergers and disturbed profiles) 
suggest that "relaxed" halos should have w 10% larger 
average c and smaller scatter. For massive systems 
(M > (few) x 10 14 M Q ), the outlook for analyzing a sam- 
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pie with a well-determined, observational selection func- 
tion is excellent, because there exist several well-defined 
catalogs o f the brightest, most massi ve clusters in X- 
rays (e.g., iReiprich fc B ohringcr 2003), of which many 
of the systems have been observed either with Chandra 
or XMM. For lower masses the number of X-ray cata- 
logs of complete samples is small, and the number with 
good coverage with Chandra or XMM observations is 
even smaller. For this purpose we have scheduled obser- 
vations of a complete, X-ray flux-limited sample of 15 
systems in the approximate mass range 10 13 — 10 14 Mq 
with Chandra. 

6.2. Hydrostatic Equilibrium 

The determination of the mass distribution from X- 
ray observations requires that hydrostatic equilibrium 
is a suitable approximation. This approximation has 
been tested for massive clusters by comparing to re- 
sults obtained from gravitational lensing (e.g., iBuotd 
120031 and references therein). There is very good agree- 
ment between the methods, especially outside the inner 
cores. This agreement is especially encouraging since 
some of the cluster s are manifestly not completely re- 
laxed (e.g., A2390, lAllen et al.ll2001l) . At the low-mass 
end, good agreement between X-ray mass measurements 
with stellar dynamics in elliptical galaxies provides fur- 
ther indication that the hydrostatic equilibrium approx- 
imation is quite accurate for obviously rel axed systems 
([Humphrey et al.ll2006l : iBridges et alJl2006t ). 

For over ten years cosmological hydrodynamical simu- 
lations have found that the hydrostatic equilibrium ap- 
proxi mation is quite accurate in massive galaxy clusters 
(e.g.. lTsai eFalfll99l iBuote fc Tsail fl^lE^rard et all 
Il996t iMathiesen et al.lll999fh Themost recent studies 
conclude that X-ray mass estimates of the most mas- 
sive, relaxed clusters should typically underestimate the 
mass by a small amount (« 10%) because of turbulent 
pressure in the hot gas, with less of an effect in lower 
mass system s (e.g.. IDolag et all 120051: iRasia et al.l [20061: 
iNagai et al.ll200lf . 

The systems analyzed in our present investigation were 
selected to be highly relaxed - as indicated by regu- 
larly shaped X-ray image morphology. This is reinforced 
by our analysis of the X-ray c — M relation (<J5|), es- 
pecially by the small intrinsic scatter in the concen- 
trations which is a robust pred i ction of CDM mod- 
els (IJhigl [20001: IDolag et~aTl 12001 IWechsler et all 12001 ; 
IMaccio et al.ll2006l) . 

Underestimating the virial masses of our relaxed clus- 
ters by 10%, as suggested by CDM simulations, would 
indicate that our measurement of C14 should be raised 
by a factor l.l -a ~ 1.02, using a = —0.17. In addition, 
since c oc (M-v^) 1 / 3 /r s , where r s is the NFW scale ra- 
dius, the total increase in C14 should be ~ 5%, provided 
the estimate of r s is unaffected by the presence of turbu- 
lent pressure. This systematic error, if real, would have 
the effect of increasing the discrepancies between the X- 
ray c — M relation and (most) CDM models. However, 
essentially all of conclusions in S}5] remain unchanged if 
we consider the bias for relaxed early forming systems 
to be 15% rather than 10% with i n the range obtained 
by current simulat ions (Jing 20001: IWechsler et "all 120021 : 
IMaccio et all l2006). The exceptions are the OCDM and 
DECDM models (g5H see Fig. gj), for which the w 5% 



boost in C14 would bring the models into slightly better 
agreement with the data. However, the OCDM model is 
still rejected at > 99.99% confidence, and the Cl m « 0.3 
open CDM models are still excluded (> 99% conf.), the 
latter provided that at least a 5% bias (of the expected 
« 10%) for relaxed, early forming systems applies. We 
conclude, therefore, that the level of systematic error sug- 
gested by CDM simulations to affect the X-ray mass mea- 
surements does not change the conclusions of our present 
study. 

6.3. Semi- Analytic Model Predictions of c — M relation 

In order to use the X-ray c — M relation for precision 
constraints on cosmological parameters, it is necessary to 
be able to predict halo concentrations produced in differ- 
ent cosmological models with high precision. This must 
be achieved with a semi-analytic procedure, because it 
is not feasible to resort to N-body simulations to fully 
investigate parameter space for obtaining confidence re- 
gions. The procedure proposed by B01 is currently one 
of the most promising models of this kind. With just 
two parameters F and K (see Sj4| it can reproduce the 
results of CDM N-body simulations with ranges of power 
spectrum shapes, as normaliza tions, matter conten t, and 
dark energy variables w (B01. iKuhlen et al.ll2005h . 

Recently, IMaccio et al.l (|2006h show that the normal- 
izations of the B01 models (i.e., K values) obtained by 
different investigators analyzing halos covering approxi- 
mately the same mass range (< 10 13 /i _1 AfQ) can differ 
by 10-2 0%. For halos with M « 10 13 - lO^/i^M©, the 
halos of IMaccio efal] (120061) pre fer F = 0.001 and K = 

2.6 while those iTasitsiomi et alj (|2004 prefer K = 3.1. 
The "cluster" halos (> 3 x 1O 14 M ) studied by D04 sug- 
gest a higher normalization K = 3.5. Because so few 
studies have investigated the high-mass halo regime, it 
is unclear whether these differences reflect numerical de- 
tails in the simulations or demand a revision of the B01 
model. More studies of the concentrations of the highest 
mass halos are very much needed. 

In t his context it is important to note that I Shaw et al.l 
(2006) have recently studied the c—M relation of massive 
halos, where most of their halos have M « 1 X 10 14 M Q . 
They obtain power-law fit parameters a = —0.12 ± 0.03 
and C14 = 8.30 ± 0.04 for the model ACDM1 but with 
(78 = 0.95. Converting their results to erg = 0.90 yields 
a 6% reduction in the concentration, C14 = 7.81 ± 0.04, 
which is just 2% larger than the D04 value quoted in 
Table [3] We consider this quite good agreement, consid- 
ering the differences between theoretical studies noted 
above for lower mass halos. 

Finally, we mention that if we adopt the lowest nor- 
malization quoted in the literature (F = 0.001, K = 2.6, 
IMaccio eT al. 2006) for our entire mass range, then in or- 
der for the ACDM1 model to match the X-ray data, it 
is necessary to boost the concentrations produced in the 
ACDM1 model by another 16% over the 10% attributed 
to the preferential selection of relaxed, early forming sys- 
tems. This 16% increase can be achieved by increasing 
either a^{~ 1.0), w (« —0.8), or both. We note that the 
99% upper limit derived in fj5j] would increase as from 

1.07 to 1.15 if K = 2.6. 

6.4. Gas Physics & Adiabatic Contraction 
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The published CDM predictions for the c — M re- 
lation we have considered in this paper are all de- 
rived from dissipationless N-body simulations contain- 
ing only dark matter. While the details of gasdynam- 
ics should be unimportant for massive systems (M > 
1O 14 M0), the effects of dissipation and feedback from 
star formation and AGN likely influence the dark mat- 
ter profile inferred from observations of lower mass ha- 
los. It is noteworthy that the systems with M ~ 
1O 13 M show the largest deviations from the ACDM1 
model ( ij5.ip . For most of these the Chandra data re- 
quire a signific ant contribution of stellar mass from the 
centr al galaxy (Humph rey et alJl2006t iGastaldello et al.1 
20061). However, allowing for adiabatic cont raction (e.g., 



zvvb). However, allowin g tor adiabatic contraction (e.g., 
Blumenthal et alj|1986t iGnedin et aT1l2004l ) of the dark 



matter profile in most cases degrades the fits, suggesting 
a more complex in terplay between the b a ryons and the 



dark matter (e 



12004 IDutton et at 



Loeb fc Peebles! l2003t lEl-Zant et all 
I2006D . Some cosmological simula- 



tions with gas predict a small, but significant, system- 
atic i ncrease (~ 3%) in the concentration of the total 
mass (|Lin et al.l |2006[ ) , but the inability of CDM sim- 
ulations to reproduce observed X-ray temperature and 
density profiles w i thin 50-100 kpc of clust e r centers (e.g., 
Lewis et al.l 120001: iMuanwong et al.l 120021 : iBorgani et all 
20041 ) makes it difficult to interpret the reliability of such 
results. Future precision cosmology studies which aim 
to use the X-ray c — M relation for low-mass systems 
(M < 10 13 M Q ) will require better understanding of the 
influence of the central galaxy on the inferred dark mat- 
ter distribution. 

7. CONCLUSIONS 

We present the concentration (c)-virial mass (M) re- 
lation of 39 galaxy systems ranging in mass from indi- 
vidual early-type galaxies up to the most massive galaxy 
clusters, (0.06 - 20) x 10 14 M o . We selected for anal- 
ysis the most relaxed systems possessing the highest 
quality data currently available in the Chandra and 
XMM public data archives. Measure ments for 24 sys- 
tems were taken from our rece n t wor k ( Hum phrey et al.l 
120061: IGastaldello et al.l 120061 : IZappacosta et all l2006h 
which populate the lower mass portion of the sample 
(M < 10 14 Af Q ). We obtained results for 15 mas- 
sive galaxy clusters in our sam ple from the studies b y 
IPointecouteau et alj (|2005l) and IVikhlinin et all (|2006l) . 

Our principal objective is to measure the c — M rela- 
tion accurately from galaxy to cluster scales and deter- 
mine whether c decreases with increasing M as generally 
predicted by CDM models. However, we also use the 
c — M relation to provide an initial demonstration of the 
ability of the c — M relation to constrain cosmological 
parameters, which also serves to highlight key sources of 
systematic error - both in the theoretical models and the 
observations. 

We parameterize the X-ray c — M relation using a sim- 
ple power-law model. The best estimates of the param- 
eters - the slope (a) and normalization (cu) evaluated 
at M = M14 = 1O 14 /i _1 M0 - were obtained via linear 
regression in log space taking into account the uncertain- 
ti es on both c and M. We em ployed the BCES method 
of lAkritas fe B crshadv (1991) for this analysis. 

Fitting the power-law model to the entire sample yields 
a = -0.172±0.026 and c 14 = 9.0±0.4 (quoted errors la). 



The slope a is negative and inconsistent with at 6.6o\ 
The previous studies of gal axy clusters (M > lO^M^) 
with Chandra and XMM b y IPointecouteau et all (|2005l ) 
and IVikhlinin et alj (|2006f ) did not place strong con- 
straints on a, and were quite consistent with a = 0; i.e., 
it is the lower mass galaxy groups that require a < 
(see IGastaldello et al.l 12006). Recent optical studies of 
the c — M relation in the group-cluster regime (see JQ) 
also do not place strong constraints on the slope and 
are very consistent with a — 0. Our analysis, therefore, 
provides crucial evidence that c de creases with increas- 
ing M, as expected in CDM models (|Navarro et al.lll997l : 
Ijrnll2000l : B01; D04: lMaccio"eTai1l2006l ). 

We compare the X-ray data to the ACDM1 model 
(with fi m = 0.3, <t$ — 0.9, n s = 1), the "concordance 
model" effectively representing the combined constraints 
from the first year of WMAP CMB data, supernovae, 
and galaxy surveys (e.g.. ISpergel et ai1l2003l ). We judge 
the median c — M relation of the ACDM1 model to be 
consistent with the empirical power-law fit provided the 
X-ray sample consists of the most relaxed, early forming 
systems, for which a systema tic increase in the concen- 
trations of ~ 10% is expe cted ([Ting 2000; Wec hsler et al.1 
l200llMacci6 et al.ll2006h . We measure an intrinsic scat- 
ter in log 10 (1 + z)c of 0.102 ± 0.004 in excellent agree- 
ment with the prediction of CDM simulations for the 
most relaxed, early forming systems. The amount of 
scatter is a robust prediction of CDM variants and pro- 
vides additional evidence that our sample comprises the 
most relaxed systems that are best suited for X-ray mass 
measurements requiring approximate hydrostatic equi- 
librium. 

The X-ray c — M relation places interesting constraints 
on o"8- Within the context of the concordance model 
noted above, the c—M relation requires 0.76 < erg < 1.07 
(99% conf.), assuming a 10% upward bias in the concen- 
trations for early forming systems. This confidence range 
is conservative as explained in tJ5.ll 

Next wc compare the X-ray data to the ACDM3 model 
(with Q m = 0.24, h = 0.73, a s = 0.76, n s = 0.96) which 
effectively updates the ACDM1 mo del using the analysi s 
of the third year of WMAP data (jSpergel et al.ll2006|) . 
The X-ray c — M relation rejects the ACDM3 model at 
the 99.99% confidence level. The primary reason for the 
poor performance of the ACDM3 model is the low value 
of erg, but the lower value of f2 m and the tilt of the power 
spectrum also contribute to the poor fit. For this com- 
parison we have assumed a uniform bias for relaxed, early 
forming halos of « 10%, as suggested by numerical sim- 
ulations. This bias would have to be « 50% to bring the 
ACDM3 model into agreement with the X-ray data. 

While the early-type galaxy and group-cluster mass 
halos (> 1O 13 M ) studied here apparently prefer 
slightly higher concentrations than predicted for typ- 
ical halos in the concordance ACDM mo del, the op- 
posit e is true for late- type galax i es fe.g.. lAlam et ail 
20021: Ivan den Bosch et al.l 120031: IDutton et al.l 120051 : 



Kuzio de Narav et all 12006: Gn edin et al1l2006l ). This 



may indicate that a selection/formation-time bias op- 
erates across the galaxy type spectrum, with late-type 
galaxies inhabiting t he low-concentrat i on tail of the dis- 
tribu tion (though see lNapolitanoll200l : lNapolitano et al.1 

Hoop). 

Since D04 provide results of power-law fits to the c—M 
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relations of massive clusters formed in a variety of CDM 
models, including an open model and several dark energy 
models, we analyze separately the X-ray c — M relation 
for the 22 systems in our sample with M > 1O 14 M . 
In this mass range we obtain cu and a ~ values con- 
sistent with those inferred from previous Chandra and 
XMM studies in the clust e r mass regime noted above 
(|Pointecouteau et all [20051 IVikhlinin et all I2006D . As 
also found for the entire sample, the c — M relation of 
the high-mass subsample is consistent with the ACDM1 
model and very inconsistent with the ACDM3 model. 

However, an open model with fi ro » 0.3 is ruled out 
(> 99% confidence) from joint consideration of only the 
X-ray c — M relation and published constraints on erg 
from the analysis of the abundances of galaxy clusters 
(erg > 0.70, see §5.2p . Since cluster abundance studies 
also find as ~ 0.45 if il m ~ 1, using the X-ray c — M data 
we find that we can also reject CDM models with f2 m w 1 
at a high significance level (> 99.9% confidence). Conse- 
quently, the combination of the X-ray c — M relation and 
cluster abundances (and local constraints on the age of 
the universe) provides novel evidence for a flat, low-f2 m 
universe with dark energy using observations only in the 
local (z « 1) universe. 

If the values of as, O m , and n s of the ACDM3 model 
are correct, agreement with the X-ray c — M relation 
may be achieved by increasing the dark energy equa- 
tion of state parameter w. We find that a quintessence 
model with w « —0.8 and as = 0.82 performs nearly as 
well as the ACDM1 model, and the larger value of w re- 



mains marginally consistent with supernova constraints 
(|Astier et al.|[2006h . 

Finally, we discuss key sources of systematic error asso- 
ciated with both the X-ray measurements and theoretical 
models that need to be addressed before the X-ray c — M 
relation is suitable for precision cosmology. In particular, 
if the virial masses are systematically underestimated by 
~ 10%, as suggested by CDM simulations, then we esti- 
mate that C14 is increased by ~ 5%, less than the w 10% 
increase expected for selecting relaxed, early forming sys- 
tems. This level of systematic error does not change the 
conclusions of our present study (see H6.2p . but it will be 
important for future studies of precision cosmology. 
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